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The phylogenetic relationships of Entada (Fabaceae) were examined to determine if two distinct 
species occur in Japan. Nucleotide sequences of the trnK and trnL introns and the intergenic 
spacer of trnL-trnF of cpDNA resolved 19 haplotypes among 39 widely collected samples of 
E. phaseoloides-Wke plants. In the resultant tree, two distinct clades corresponded to two groups 
also distinguished by diagnostic seed characters. The clades agree with the expected chorology of 
samples and with the diagnostic morphology from previous taxonomic work. One clade corresponds 
to Entada phaseoloides and the other to E. tonkinensis. A distinctive concentration of haplotypes 
was observed. Although seeds of Entada phaseoloides are renowned as sea drift seeds, our results do 
not suggest frequent gene flow caused by long distance seed dispersal. 
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Entada Adans. (Leguminosae: Mimosoi- 
deae) consists of about 30 species (Nielsen 1981, 
1992, ILDIS 2005) distributed in the tropics and 
subtropics. While more than half of the species 
in Entada are shrubs or subshrubs with small 
pods, some species are lianas with gigantic pods 
that can be more than 1 m long. The large woody 
seeds of Entada are also known for drifting in sea¬ 
water. These remarkable seeds are called “match¬ 
box beans” (the common name for Entada phase¬ 
oloides (L.) Merr. and E. rheedii Spreng.) or “sea 
hearts” ( E. gigas (L.) Fawc. & Rendle), and are 
frequently found on beaches far from their origin. 

Entada phaseoloides is widely distributed in 


tropical and subtropical areas of the Pacific, from 
eastern to southeastern continental Asia, the south 
Asian islands, Oceania, and Polynesia. The north¬ 
eastern limit of the species’ range is Japan (Nielsen 
1992, Huang & Ohashi 1993). The wide distribu¬ 
tion of Entada phaseoloides is thought to have 
been achieved by long distance dispersal of the 
seeds. Although Entada phaseoloides has some 
distinctive features in its fruits and seeds, consid¬ 
erable taxonomic confusion exists as to its clas¬ 
sification in Japan and adjacent areas. While most 
studies treat the species of Entada in Japan as E. 
phaseoloides, some works have variously treated 
it as E. rheedii Spreng., E. parvifolia Merr., or E. 
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koshunensis Hayata & Kaneh. Some authors rec¬ 
ognize two distinct species (Ho 1985, Huang & 
Ohashi 1993, Ohashi 2001), but the recognition of 
E. rheedii and E. parvifolia in the Ryukyu Islands 
(Okinawa), Japan, is apparently misleading and E. 
koshunensis cannot readily be separated from E. 
phaseoloides according to a detailed taxonomic 
study that focused on this taxonomic confusion 
(Tateishi et al. 2008). 

Recently, Wakita et al. (2005) surveyed the 
number of leaflet pairs on the second pinnae, 
morphological characters of the seeds, and single 
strand conformation polymorphism (SSCP) band 
patterns of the rpsl6 intron in chloroplast DNA 
(cpDNA) for species of Entada in the Ryukyu 
Islands. Although the number of pairs of leaflet 
on the second pinnae did not clearly divide the 
samples into two groups, the seeds and molecu¬ 
lar markers clearly indicated the two groups of 
Entada. The first group has a small (‘S’), convex 
seed type, 29.4^16.9 mm wide by 32.4-52.7 mm 
long with an angular margin. The second group 
has a large (‘L’), compressed seed type, 39.5-59.1 
mm wide by 45.8-64.7 mm long with a rounded 
margin (Wakita et al. 2005). Different SSCP band 
patterns, caused by a nucleotide substitution and 
some indels in the rpsl6 intron, were also clearly 
observed in the two groups. Wakita et al. (2005) 
considered the first group as corresponding Enta¬ 
da koshunensis and the second to E. phaseoloides 
based on seed characteristics, following the cri¬ 
teria of Ho (1985) for distinguishing species of 
Entada in Taiwan. Wakita et al. (2005), however, 
did not resolve the taxonomic status of these two 
groups due to insufficient sampling limited only 
to the Ryukyu Islands. A molecular phylogenetic 
study using additional samples from a wider geo¬ 
graphic area is necessary to clarify the presence of 
two species in Japan. 

In this paper, we report the results of molecu¬ 
lar phylogenetic analyses using samples obtained 
from a broad geographical range, including the 
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presumed type locality of Entada phaseoloides 
in Indonesia. Key seed shape characters used in 
Wakita et al. (2005) were also examined to deter¬ 
mine the correspondence between seed character 
groups and molecular phylogenetic groups. The 
goal of this study was to reveal the phylogenetic 
relationships of Entada phaseoloides and closely 
related species, and to determine if indeed two 
species of Entada occur in Japan. 

Materials and Methods 

Entada phaseoloides-\ike plants were col¬ 
lected from eastern to southeastern Asia and the 
South Pacific (Table 1). Samples from Indone¬ 
sia and Vanuatu were obtained from individu¬ 
als cultivated in the Bogor Botanical Garden 
(Bogor, Indonesia) and the Tsukuba Botanical 
Garden (Ibaraki, Japan), respectively (Table 1). 
Since plants of Entada are generally very large, 
spreading over the forest canopy, distinguishing 
individuals was difficult. We collected more than 
one sample per population only when individu¬ 
als could be clearly distinguished. Seed samples 
were not obtained from all the populations. Key 
seed shape characters of all available seeds were 
observed, and seed type (S or L) was determined 
according to Wakita et al. (2005). 

Leaf samples of the species of Entada 
obtained from eastern and southeastern Asia, 
Oceania, and the South Pacific were used for a 
molecular phylogenetic study (Table 1). Entada 
gigas (L.) Fawc. & Rendle, E. rheedii , and E. 
spiralis Ridley, which belong to the same subsec¬ 
tion ( Entada subsect. Entada) as E. phaseoloides, 
were also included in the phylogenetic analysis. 
Outgroups were Entada abyssinica Steud. (from 
Tanzania) and E. polystachya DC. (from Mexico). 
Those species were confirmed as not being in the 
ingroup by our preliminary analysis (N. Wakita et 
al., unpublished data) using published trnL-trnF 
sequences of Mimosoideae (Luckow et al. 2000, 
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2003). A total of 47 samples, including outgroups, 
was used for the sequence analysis. Voucher 
specimens of the samples used in this study were 
deposited in QCNE, TI, TNS, or URO (Table 1). 

Prior to genomic DNA extraction, crushed 
leaf tissue was washed using HEPES buffer fol¬ 
lowing the methods of Setoguchi & Ohba (1995). 
Total genomic DNA was extracted from washed 
leaf pellets using the methods of Doyle & Doyle 
(1987). The two regions of cpDNA, the trnL 
intron including the trnL-trnF intergenic spacer 
(IGS; hereafter referred to as trnL-trnF) and the 
trnK intron (including matK\ hereafter referred to 
as trnK), were amplified by PCR. The following 
oligonucleotide primers were used for PCR ampli¬ 
fication and sequencing: c, d, e, and f for trnL- 
trnF, designed by Taberlet et al. (1991); trnK- F 
(forward: 5 -TGGGTTGCTAACTCAACGG-3') 
and trnK- R (reverse: 5 -GGAACTAGTCGGAT- 
GGAGT-3') designed according to the complete 
cpDNA sequence data of tobacco (GenBank 
accession numbers Z00044 and S54304), matK-10 
(Kato et al. 1998), matK-AF (Takayama et al. 
2005), matK-1777L, and matK-1932R (Hu et al. 
2000) for trnK. The PCR reaction mixture con¬ 
tained 2.0 pL of 0.2 mmol/L dNTP solution, 1.0 
pL of 5 mmol/L of each primer, 0.5 unit of ExTaq 
DNA polymerase (TaKaRa Bio, Ohtsu, Shiga, 
Japan), 2.5 pL of 10x ExTaq buffer, and 10-30 ng 
of genomic DNA in a total volume of 25 pL. PCR 
amplification was performed in a PCR Thermal 
cycler Dice (TaKaRa Bio). The PCR cycle for 
trnL-trnF was 95°C for 1 min; followed by 25 
cycles at 95°C for 45 s, 51°C for 1.5 min, and 72° 
C for 1.5 min; and a final extension at 72°C for 10 
min. For trnK, the settings were 96°C for 1 min; 
followed by 35 cycles at 96°C for 45 s, 48°C for 
1 min, and 72°C for 1.5 min; and a final extension 
at 72°C for 15 min. The extension time in both 
PCR programs was increased by 2 s every cycle. 
Amplified DNA was purified using a GeneClean 
III DNA purification Kit (BIO 101, Carlsbad, 
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CA, USA). Purified DNA was directly sequenced 
using an ABI PRISM Big Dye Terminator ver¬ 
sion 3.1 Cycle Sequencing Ready Reaction kit 
and an ABI PRISM 377 DNA sequencer (Applied 
Biosystems, Foster City, CA, USA). Multiple 
sequence alignment was performed with ClustalX 
ver. 1.83 (Thompson et al. 1994) using the default 
setting, and the obtained alignment was manually 
modified using Se-Al Sequence Alignment Editor 
version 2.0all (Rambaut 2002). Combined nucle¬ 
otide sequences were used for further analyses 
because the incongruence length difference (ILD) 
test implemented in PAUP* ver4.0. detected no 
significant incongruence between the two genetic 
regions for phylogenetic analyses. 

Maximum-parsimony (MP) and neighbor¬ 
joining (NJ) analyses were conducted using 
PAUP* version 4.0 beta 10 (Swofford 2002). 
Branch support was assessed with 1,000 rep¬ 
licates for MP and 10,000 for NJ. A Bayesian 
phylogenetic analysis was performed in MrBayes 
version 3.1.2 (Ronquist & Huelsenbeck 2003) 
with the best-fit model (GTR+I+G) for both of 
trnL-trnF and trnK selected by Akaike’s infor¬ 
mation criterion in MrModeltest 2.2 (Nylander 
et al. 2004). Metropolis-coupled Markov chain 
Monte Carlo sampling was performed with one 
cold chain and three heated chains that were run 
for 1,000,000 generations. Trees were sampled 
every 100 generations. In total, 10,000 trees were 
obtained and the first 2,500 samples from each 
run were discarded as bum-in to ensure the chains 
had become stationary. A majority rule consensus 
tree was generated to compute posterior probabili¬ 
ties. Insertion/deletion events (indels) were not 
included in the NJ and Bayesian analyses. Indels 
of more than two base pairs were treated as binary 
(0/1) characters for the MP analysis if not caused 
by a simple sequence repeat. 


WAKITA & AL.: Entada in Japan by cpDNA Phylogeny 


Nil-Electronic Library Service 



The Japanese Society for Plant Systematics 


Acta Phytotax. Geobot. 


Vol. 59 


VO <N —• cn 

vo r-~ r- r~~- oo 

OO OO OO 00 00 

N (N G) M N 

m cn m co m 

P P P 13 D 

W W W W W 


h- iTl OO 

r- r- r- 
t ^ t 
GO 00 00 
N N M 
m co co 

P P P 

WWW 


m On oo O Ov 

O O O -H ^ 

OO 00 00 00 oo 

O) M (N (N n 

co m m co m 

P p P P P 

w w w w w 


Tf tJ- 

OO 00 00 
(N N G) 

ro co co 

P D P 

WWW 


c/5 i/i c/5 w W c/3 


W W C /3 


^rH^UO^TfpHNOOhON-H-H-th-^'-HGMfh-l'^GG^r'O'Oh 
OOOOOOOOOOOOOOOOOOOO-^ — 

ta.n.fcfcfcfcdfca.tt.HHSHtt.HHhhHA.a.fcfcfcfl.a.dA.fcfl. 


“3 =8 
^ ^ ^ 

Si S; S; 


<N *n 
S <N fN 
Os On 
Q O 

°s ^ ^ 
^ ^ 
^ S; fe; 


<N *0 w 
. <N °0 

« ^ °? 

<3 •'*» *0 
S ^ ^ 
S «N <N 
§ ^ ^ 

s Si si 


^ -S z 

J3 t 0 1 Z 


^ ^ Os 
°8 °0 »N 


*C ^ | 
w w ^ 


q q q 

O O CJ 

"O "O "O 
cj e3 ctf 

o o o 


° 3> g 

5c® 


6 § 

O cv . 

a ^ g 

JS g & 

. - i* £ 

! ■§ 

I g I 


o o o 
CO CO CQ 
u u u 
O O O 
W) W) 00 
o o o 
CQ CO CO 
o a> <u 
-c -c W 


'O T3 'TO 

<D <D <D 

# > .> .> 

3 3 3 
u u u 


.5 .S ts 
•a o Z 

H H ni 


3 3 a 

U O Q 


t: c c t t; c 

ra ra ra la an c <3 03 

p, p p n q v p ( Qh 


6 6 6 S 3 £ 6 e — — a | | 

5 3 S i § H £ > > 'S o E § !{ 


(Li c /2 «J 

^ s § S 


3 3 £h Vh q t-G c“ ?? .,_, _*2 H GO 0 > 

5oOOOOOC>£;3' , 35p3 

t%&ZZXZZ<ggi*t££oO l 


« 5 5 
■a o o 
sg o. a 
pp >—) >-J 


Nil-Electronic Library Service 




AUSTRALIA Queensland York Cape Peninsula Y. Tateishi & Omine 53109 ( URO) R03 EU328423 EU328486 

INDONESIA Java Isl. Cultivated in the Bogor Botanical Garden N.W 250510-G27 (U) R04 EU328424 EU328487 


The Japanese Society for Plant Systematics 


November 2008 


WAKITA & AL.: Entada in Japan by cpDNA Phylogeny 


187 



Nil-Electronic Library Service 



The Japanese Society for Plant Systematics 


188 Acta Phytotax. Geobot. Vol. 59 



Fig. 1. Geographic distribution of Entada haplotypes. Haplotypes are designated by name as they appear in Table 1 and Fig. 2. 
Circles and triangles indicate the clade to which each haplotype belongs (see Fig. 2). Numbers in parentheses after haplotype 
names indicate numbers of individuals. Seed type is shown as S or L according to Wakita et al. (2005). (a) -(c): Distribution 
maps, (d): Pictures showing seed types. Scale bar = 1cm. 
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Fig. 2. Strict consensus tree of 123 most parsimonious trees of E. phaseoloides-Mke samples and their allies obtained from trnL- 
IrnF and trnK sequences (length = 202 steps, Cl = 0.861, R1 = 0.8607, RC = 0.741). Measures of support are given at each 
node: NJ bootstrap support/MP bootstrap support/Bayesian posterior probability. 
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Results 

Seed morphology 

Using the key character of seed shape report¬ 
ed in Wakita et al. (2005), the seeds were sorted 
into one of two groups. The seed type assigned 
to all available seeds is shown in Table 1. Type 
S seeds were obtained from Japan (Okinawa), 
southern Taiwan, Indonesia, Australia, Vanuatu, 
and Samoa. Type L seeds were obtained from 
Japan (Kagoshima), northern Taiwan, China, and 
Vietnam. The geographic distribution of seed 
types is mapped in Fig. 1. 

Molecular analyses 

The trnL-trnF sequence was 1,050-1,057 
bp long and the trnK sequence was 2,513-2,537 
bp long. Using the combined sequences, 19 
haplotypes were resolved among the 39 Entada 
phaseoloides ASks, samples. Twenty-four indels 
were obtained across the matrix of all haplotypes 
and three were coded as potentially phylogeneti- 
cally informative for parsimony analyses. The 
phylogenetic analysis using Entada abyssinica 
and E. polystachya as outgroups produced 164 
most parsimonious trees of 226 steps, with a con¬ 
sistency index of 0.876, retention index of 0.877, 
and rescaled consistency index of 0.769 (Fig. 2). 
The NJ and Bayesian tree topologies were nearly 
consistent with the MP tree presented here (Fig. 2). 
Entada gigas was sister to another clade, which 
was further divided into two major clades. One 
of those clades was composed of Entada rheedii 
and E. spiralis (branch support (BS) by bootstrap 
or posterior probability of 71 (MP analysis), 85 
(NJ analysis), and 0.90 (Bayesian analysis); BS is 
presented in this order in the following text). The 
second clade was composed of all Entada pha- 
seoloides-like samples (BS of 96, 96, and 1.00). 
This latter clade was further subdivided into two 
clades (Clades I and II), both of which were very 


well supported (BS of 99, 99, and 1.00 for Clade 
I and 98, 99, and 1.00 for Clade II). Clade I con¬ 
sisted of samples from southern Taiwan, Vanuatu, 
Samoa, Australia, and Indonesia, where S-type 
seeds were obtained. Clade II consisted of sam¬ 
ples from northern Taiwan, Hong Kong, and Viet¬ 
nam, where L-type seeds were obtained. Samples 
collected from the type locality of Entada phase- 
oloides (Ambon Island, Indonesia) were included 
in Clade I. 

Discussion 

Presence of two species of Entada in Japan 

The two well-defined clades (I and II) in 
our phylogenetic tree are clearly consistent with 
the two seed type groups (S or L) recognized by 
Wakita et al. (2005). Our results clearly indicate 
the presence of two species of Entada in Japan. 

Entada phaseoloides was originally described 
by Merrill (1917). The type specimen (drawn as 
Faba marina in tab 4, Herbarium Amboinense 5: 
5-8, Rumphius 1747) was obtained from Ambon, 
Indonesia. According to some floristic works that 
treated this species (Walker 1976, Nielsen 1992), 
the diagnostic characters of Entada phaseoloides 
are 1-3 (or 4) pairs of asymmetric elliptic leaflets 
and seeds 3-6 cm long. Entada phaseoloides is 
distributed throughout subtropical to tropical areas 
from Southeast Asia to the Pacific islands. This 
general characterization is applicable to all sam¬ 
ples considered to be Entada phaseoloides-like in 
this study. The wide distribution reported thus far 
may represent the sum of the two species revealed 
here. Our study, however, clearly suggests distinc¬ 
tive distribution areas for the two species. Sam¬ 
ples included in Clade I were from Southeast Asia 
to the South Pacific, whereas samples in Clade II 
was collected only in areas extending from Viet¬ 
nam to Japan (Fig. 1). Since the sample obtained 
from the type locality of Entada phaseoloides 
(Ambon Island, Indonesia) appears in Clade I, 
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we conclude that this clade corresponds to E. 
phaseoloides. All available seeds obtained from 
samples in Clade I were of the S type. Thus, the 
diagnostic seed character of Entada phaseoloides 
is a convex seed with an angular margin, and the 
range of E. phaseoloides encompasses the islands 
of eastern and southeastern Asia, a part of Oce¬ 
ania, and Polynesia (Fig. 1). 

The second clade (Clade II) was composed 
of samples from Japan (Amami-Oshima and Yaku 
islands), northern Taiwan, Hong Kong, and Viet¬ 
nam. All available seeds from these samples were 
type L (Table 1, Fig. 1). The samples compris¬ 
ing Clade II have generally been misidentified in 
Japan and Taiwan as Entada phaseoloides. After 
surveying various taxonomic works from Japan, 
Taiwan, China, and Vietnam, the seed morphol¬ 
ogy of Entada tonkinensis Gagnep. appears to be 
consistent with the samples in Clade II. Entada 
tonkinensis was originally reported from Vietnam 
(Gagnepain 1911), and its characteristically large 
seeds were recorded as 50 mm wide and 60 mm 
long. Clade II contained samples from Vietnam 
and the seeds obtained from this clade were 39.5 
-59.1 mm wide and 45.8-64.7 mm long, which 
is consistent with the seed size reported for E. 
tonkinensis. Although the species has not been 
reported outside Vietnam, our study suggests a 
wider distribution extending to China, Taiwan, 
and Japan. 

Unlike the above discussion, Wakita et al. 
(2005) assigned Entada phaseoloides to the group 
with type L seeds and E. koshunensis to the group 
with type S seeds, based on samples obtained 
from the Ryukyu Islands. Their assignment of 
Entada phaseoloides to the type L seed group is 
not likely correct, as discussed above. Regarding 
their assignment of Entada koshunensis to the 
type S seed group, our phylogenetic tree does not 
show clear evidence to distinguish this species. 
Entada koshunensis was originally described by 
Hayata (1921) based on a sample collected in 
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the Henchun Peninsula, southern Taiwan. Some 
studies distinguish this species from Entada 
phaseoloides based on its smaller seeds (ca. 3 
cm) and greater number of leaflet pairs (Hayata 
1921, Ho 1985, Huang & Ohashi 1993), but other 
studies treat it as E. phaseoloides because of its 
morphological resemblance (Walker 1976, Huang 
& Ohashi 1977). In this study, four samples col¬ 
lected from the type locality of Entada koshun¬ 
ensis (Henchun Peninsula, southern Taiwan) had 
haplotypes P01 and P03. These haplotypes were 
included in Clade I and formed a clade having 
weak support with samples from the Ryukyu 
Islands (Fig. 2), but no clear morphological or 
biogeographical association with this clade was 
observed. Further taxonomic study also supports 
the difficulty in distinguishing Entada koshunen¬ 
sis as a separate species (Tateishi et al. 2008). 

Geographic distribution ofE. phaseoloides and E. 
tonkinensis haplotypes 

This study revealed the approximate distri¬ 
butions of haplotypes of Entada phaseoloides and 
E. tonkinensis (Fig. 1). Both species have very 
large seeds that can drift in seawater, which might 
contribute to their population structure. Other sea- 
dispersed plants have recently been studied using 
molecular markers and evidence of long distance 
dispersal in haplotype structure has been reported 
(Takayama et al. 2006, Nettle & Dodd 2007). 
For example, Hibiscus tiliaceus L. (Malvaceae) 
has a wide distribution throughout the Old World 
tropics. Common haplotypes were distributed 
throughout the Pacific and Indian ocean areas, and 
genetic differentiation among populations was 
very low (Takayama et al. 2006). However, our 
results from the study of Entada in Asia and the 
Pacific region are somewhat contradictory. Based 
on the distribution of haplotypes (Fig. 1), few 
widespread haplotypes in Entada phaseoloides 
are observed and a geographic concentration of 
haplotypes appears to exist. Haplotypes P01- 
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P09 occur in Japan and Taiwan, haplotypes P11— 
P14 in Indonesia, and haplotypes P10 and PI 5- 
P17 in Australia and Samoa (Fig. 1). Although the 
number of samples surveyed was limited, no com¬ 
mon haplotypes were found among these regions. 
The only case of long distance distribution of a 
common haplotype is in Australia and Samoa, 
separated by about 5,300 km, for haplotype PI7. 

The seeds of Entada phaseoloides are viable 
for at least 1 year in seawater (Nakanishi 1994); 
thus, seeds can be dispersed by ocean currents 
over very long distances. The distribution of each 
haplotype, however, was geographically limited 
(Fig. 1). One possible explanation for this is the 
low ability of seeds to establish in new localities 
after successful dispersal. Entada phaseoloides is 
rarely seen in littoral areas and is more common 
in the canopy of inland forest. Although seeds 
frequently reach open sandy beaches, subsequent 
colonization may be difficult. Other factors such 
as the direction of ocean currents may also result 
in the geographically limited haplotype distribu¬ 
tion. Further studies using more samples from 
various localities are necessary to clarify this 
issue. 

For Entada tonkinensis, Clade II included 
only two haplotypes, whereas Clade I of E. pha¬ 
seoloides contained 17 haplotypes from 23 sam¬ 
ples. This perhaps reflects the population size and 
distribution ranges of these two species. In Clade 
II, a common haplotype (T01) was found in eight 
populations in five regions (Fig. 2). The distribu¬ 
tion of haplotype T01 might have been due to dis¬ 
persal of seeds from continental Asia northward 
to Taiwan and Japan along the island chain, as the 
habitat of Entada tonkinensis in Japan and Hong 
Kong is on slopes exposed to the sea at 20-150 
m. The occurrence in Taiwan, however, is inland 
about 15-30 km from the seashore. How this spe¬ 
cies has successfully colonized an inland forest 
from the seashore on which seeds washed ashore 
is still unknown. As the vines of this species grow 
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very large, individual growth and seed production 
over many years might be one explanation for the 
colonization of an inland habitat. 

In this study, a phylogenetic analysis using 
cpDNA sequences revealed that two distinct 
phylogenetic groups of Entada are present in 
Japan. Morphologically, one clade corresponds 
to Entada phaseoloides and the other to E. tonk¬ 
inensis. Although seeds of Entada phaseoloides 
are remarkable examples of sea-drifted seeds, evi¬ 
dence of frequent long distance seed dispersal was 
not obtained. Even if the seeds can be dispersed 
by ocean currents across very long distances 
(Ridley 1930), the geographical concentration of 
haplotypes suggests that another factor, such as 
difficulty in establishment, prevents greater gene 
flow among regions. 
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